function model = ldaclassify(stego, cover)
    [ns, nc, space] = LDA(stego, cover);
    model = [];
    model.space = space;
    model.stego_mean = mean(ns);
    model.cover_mean = mean(nc);
    % ab = 0 : stego < threshold  cover > threshold
    % ab > 0 : stego > threshold  cover < threshold
    model.ab = (model.stego_mean > model.cover_mean);
    model.threshold = mean([ns;nc]);
    %model.threshold = (model.stego_mean + model.cover_mean)/2;
end
